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1. Introduction 

One of the major goals at the LHC is the production of Higgs boson(s) [1]. In the Minimal 
Supersymmetric Extension of the Standard Model (MSSM) two complex Higgs doublets are intro- 
duced to give masses to up- and down-type fermions [2]. After electroweak symmetry breaking 
there are five physcial Higgs states, two CP-even neutral Higgs bosons h^H, one neutral CP-odd 
Higgs state A and two charged Higgs bosons H^. At tree level, the Higgs sector can be param- 
eterized by two independent parameters, the pseudoscalar Higgs boson mass Ma and the ratio of 
the two vacuum expectation values (VEV) of the two complex Higgs doublets, tanjS = V2/V1. The 
Higgs couplings to quarks and gauge bosons are modified with sin and cos of the mixing angles a 
and j8 with respect to the Standard Model (SM) couplings, where a denotes the h,H mixing angle. 
The bottom (top) Yukawa couplings are enhanced (suppressed) for large values of tan j3 , so that top 
Yukawa couplings play a dominant role at small and moderate values of tanjS. 

At the LHC and Tevatron neutral Higgs bosons are copiously produced via gluon fusion gg — )• 
h,H,A, which is mediated in the case of h,H by (s)top and (s)bottom loops [3]. The pure QCD 
corrections to the (s)quark loops have been obtained including the full Higgs and (s)quark mass 
dependences and increase the cross sections by ~ 100% [4]. This result can be approximated 
by very heavy top (s)quarks with ~ 20 — 30% accuracy for tanj3 < 5 [5]. In this limit the next-to- 
leading order (NLO) QCD [6] and later the next-to-next-to-leading order (NNLO) QCD corrections 
[7] have been obtained, the latter leading to a moderate increase of 20-30%. Finite top mass 
effects at NNLO have been discussed in [8]. Finally, the estimate of the next-to-next-to-next-to- 
leading order effects [9] indicates improved perturbative convergence. The full supersymmetric 
(SUSY) QCD corrections have been obtained in the limit of heavy SUSY particle masses [10] and 
more recently including the full mass dependence [11]. The electroweak loop effects have been 
calculated in [12]. In this article we will describe in Section § the calculation of the full SUSY- 
QCD corrections in gluon fusion to h,H, and we will present for the first time numerical results 
for the total cross section. In Section |3| we will discuss the consistent derivation of the effective 
Lagrangian for the scalar Higgs couplings to gluons after the gluino decoupling. 

2. Gluon Fusion 

At leading order (LO) the gluon fusion processes gg ^ h /H are mediated by heavy quark and 
squark triangle loops, cf. Figjl], the latter contributing significantly for squark masses < 400 GeV. 
The LO cross section in the narrow- width approximation can be obtained from the h/H gluonic 
decay widths, [3,13] 
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Figure 1: Diagrams contributing to —>/;,// at leading order. 



where T;,/^ = Mf^j^/s with s being the squared hadronic cm. energy and Tg^g = Am^^^^/Mj^^^. 
The LO form factors are given by 
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And the gluon luminosity at the factorization scale Hf is defined as 

ax Jz X 

where g{x,pLp) denotes the gluon parton density of the proton. The NLO SUSY-QCD corrections 
consist of the virtual two-loop corrections, cf. FigJH, and the real corrections due to the radiation 
processes gg — > gh/H,gq — > qh/H and qq —?■ gh/H, cf. Fig.^. The final result for the total hadronic 
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Figure 2: Some generic diagrams for the virtual NLO SUSY-QCD corrections to the squark contributions 
to the gluonic Higgs couplings. 



cross sections can be split accordingly into five parts, 



G{pp^h/H+X) = o^ 
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The strong coupling constant is renormalized in the MS scheme, with the top quark and squark 
contributions decoupled from the scale dependence. The quark and squark masses are renormalized 
on-shell. The parton densities are defined in the MS scheme with five active flavors, i.e. the top 
quark and the squarks are not included in the factorization scale dependence. After renormalization 
we are left with coUinear divergences in the sum of the virtual and real corrections which are 
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absorbed in the renormalization of the parton density functions, so that the result Eq. ( |2.4| ) is finite 
and depends on the renormalization and factorization scales /Xr and jip, respectively. The natural 
scale choices turn out to be = Hf Mj,/fj- The numerical results are presented for the modified 




Figure 3: Typical diagrams for the real NLO QCD corrections to the squark contributions to the gluon 
fusion processes. 



small aeff scenario [14], defined by the following choices of MSSM parameters [m, = 172.6 GeV], 
Mg = 800GeV tan/3 =30 

Mg = 1000 GeV /i =2 TeV (2.5) 

M2 = 500 GeV Ai=A, = -1.133 TeV. 

In this scenario the squark masses amount to 

nif- = 679 GeV m^^ = 935 GeV 
= 601 GeV m-^^ = 961 GeV . 

Fig. displays the genuine SUSY QCD corrections normahzed to the LO bottom quark form factor. 
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Figure 4: The genuine SUSY QCD corrections normalized to the LO bottom quark form factor. Real 
corrections: red (light gray), virtual corrections: blue (dark gray), compared to the A/, approximation (dashed 
lines). Ah has been renormahzed in the MS scheme. 



i.e. Ay {ih) — ('^ft)(l +^sf/5y^)- The corrections can be sizeable, but can be described 
reasonably with the usual At, approximation [15], if At, is renormalized in the MS scheme. 

3. Decoupling of the Gluinos 

In this section we will address the limit of heavy quark, squark and gluino masses, where 
in addition the gluinos are much heavier than the quarks and squarks. For the derivation of the 
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effective Lagrangian for the scalar Higgs couplings to gluons we analyze the relation between the 
quark Yukawa coupling Xq and the Higgs coupling to squarks Xq in the limit of large gluino masses. 
We define these couplings at leading order in the case of vanishing mixing, 



Art = 



with K = 2 
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(3.1) 



where g-Q denotes the normalization factor of the MSSM Higgs couplings to quark pairs with re- 
spect to the SM. In the following we will sketch how the modified relation between these couplings 
for scales below the gluino mass Mg is derived. For details, see Ref. [16]. We start with the unbro- 



ken relation between the running MS couplings of Eq. ( |3.1[ ) and the corresponding renormalization 
group equations (RGE) for scales above Mg. If the scales decrease below Mg the gluino decouples 
from the RGBs leading to modified RGBs which are different for the two couplings Xq and kXq 
so that the two couplings deviate for scales below Mg. The proper matching at the gluino mass 
scale yields a finite threshold contribution for the evolution from the gluino mass scale to smaller 
scales, while the logarithmic structure of the matching relation is given by the solution of the RGBs 
below Mg. In order to decouple consistently the gluino from the RGB for gluino mass scales large 
compared to the chosen renormalization scale, a momentum substraction of the gluino contribution 
for vanishing momentum transfer has to be performed [17]. We refer the reader to [16] for details 
and give here directly the result for the modified relation between the quark Yukawa coupling and 
the effective Higgs coupling to squarks taking into account the proper gluino decoupling: 
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where iuq is the pole mass and MO denotes the momentum substracted coupling, which is taken at 
the squark mass scale, which is the proper scale choice of the effective Higgs coupling to squarks 
and which is relevant for an additional large gap between the quark and squark masses. 

Taking into account the radiative corrections to the relation between the effective couplings af- 
ter decoupling the gluinos leads to the following effective Lagrangian in the limit of heavy squarks 
and quarks. 
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^■^e,Mo('"e)/'"|- The cofficient Csqcd is given by 

^SQCD — -T ■ 



(3.4) 



It is well-defined in the limit of large gluino masses and thus fulfills the constraint of the Appelquist- 
Carazzone decoupling theorem [18]. 

4. Conclusions 



We have presented first results for the NLO SUSY QCD corrections to gluon fusion into 
CP-even MSSM Higgs bosons, including the full mass dependence of the loop particles. The 
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genuine SUSY-QCD corrections can be sizeable. We furthermore demonstrated, that the gluino 
contributions can be decoupled in the large M| Umit in accordance with the Appelquist-Carazzone 
theorem. 
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